A new structural model for the Si(331)-(12xl) reconstruction 
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A new structural model for the Si(331)-(12x 1) reconstruction is proposed. Based on scanning tun- 
neling microscopy images of unprecedented resolution, low-energy electron diffraction data, and first- 
principles total-energy calculations, we demonstrate that the reconstructed Si(331) surface shares 
the same elementary building blocks as the Si(110)-(16x2) surface, establishing the pentamer as a 
universal building block for complex silicon surface reconstructions. 



The study of semiconductor surface reconstructions 
has been an area of active research for many years 
and has gained tremendous importance with the advent 
of low- dimensional heteroepitaxial semiconductor nano- 
structures such as quantum dots and quantum wires []J . 
The creation of a surface results in broken bonds, called 
dangling bonds. Dangling bonds are energetically un- 
favorable causing surface atoms to rearrange or recon- 
struct. This often results in highly complex atomic struc- 
tures, whose determination remains a formidable chal- 
lenge and requires the complementary role of different 
experimental and theoretical methods. 
In order to lower the surface energy, silicon surfaces adopt 
a variety of strategies allowing to reduce the number of 
dangling bonds. Despite the large number of known sur- 
face reconstructions, one frequently encounters common 
elementary structural building blocks [2, 3]. Identifying 
these building blocks is important not only for a better 
understanding of these surfaces, but also serves as a guide 
for the elaboration of new structural models. 
Two of the most important strategies, encountered for 
instance on Si(lOO) 3 ^^^ Si(lll) Q, are respectively 
the formation of dimers, where two surface atoms pair up 
to eliminate their dangling bonds, and the appearance of 
adatoms, which bond to three surface atoms thus satu- 
rating three dangling bonds. An important step towards 
the understanding of high-index group IV surfaces with a 
surface normal in between the (111) and (100) direction 
was the introduction of an additional reconstruction ele- 
ment by Dabrowski et al. [6]. They proposed a six- fold 
coordinated surface self-interstitial which is captured by 
a conglomerate of surface atoms [3, [8| . This concept was 
subsequently adapted by An Q and theoretically ana- 
lyzed by Stekolnikov [lO, HH to explain the pairs of pen- 
tagons observed in scanning tunneling microscopy (STM) 
images of the reconstructed Si(llO) surface. 

In this letter we focus on the atomic structure of 
the Si(331)-(12xl) reconstruction. We present high- 
resolution STM images resolving for the first time rows of 
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FIG. 1: (Color online) a) Side cutaway showing the crystal lat- 
tice of silicon in the (110) plane. The dashed line follows the 
bulk-terminated surface for several important orientations, b) 
Experimental LEED pattern at 35 eV beam energy. The po- 
sitions of missing spots are indicated by the white circles, c) 
Sketch of the LEED pattern with (1x1) (red) and (12x1) 
(blue) reciprocal unit cells and spot labels (black). The bulk 
directions are also given. The positions of the missing spots 
are indicated by empty black circles. The orientation of the 
gUde plane is indicated by a dashed line. 



pentagons very similar to the ones observed on Si(llO)- 
(16x2). Si(331), whose surface normal is located 22.0^ 
away from the (111) direction towards (110) (see Fig. 
d^)), is an important surface, since it is the only con- 
firmed planar silicon surface with a stable reconstruc- 
tion located between (111) and (110). Since the discov- 
ery of the Si(331)-(12xl) reconstruction more than 17 
years ago [12] several structural models containing dimers 
and adatoms have been proposed [l^, [ij]. However, 
none of these models is able to explain the pentagons 
observed in our STM images. Combining the comple- 
mentary strength of STM, low-energy electron diffraction 
(LEED), and first-principles total-energy calculations, we 
derive a new structural model containing surface self- 



interstitials as basic building blocks. 
Sample preparation and experiments were carried out in 
a ultra-high vacuum chamber with a residual gas pres- 
sure below 3 X 10~^^ mbar equipped with an Omicron 
LT-STM and Omicron Spectaleed LEED/Auger optics. 
Boron doped Si(331) samples from CrysTec with a re- 
sistivity of 0.1-30 Q cm were slowly degassed by direct 
current heating up to 1060^C while keeping the pressure 
below 5 X 10~^^ mbar, followed by subsequent cleaning 
via repeated flashing to 1260^0 and slow cooling across 
the (1x1) to (12x1) phase transition at 810^C [l2|. This 
procedure gives an almost perfectly ordered surface with 
a very low number of defects. Special care was taken 
in order to avoid Ni contamination. STM measurements 
were performed at 77 K using etched W tips. 
The theoretical results are based on first-principles cal- 
culations performed with ABINIT [l5[. This code em- 
ploys a density-functional theory scheme within a plane- 
wave expansion, where we have used norm-conserving 
pseudopotentials and the local-density approximation for 
the exchange-correlation potential. ABINIT has already 
been successfully applied to describe the Si(113) surface 
[16]. The present (331) surface has been studied using a 
slab super cell containing 240 bulk silicon atoms, plus a 
variable number of surface atoms (depending on the re- 
construction). The bottom surface was passivated with 
36 hydrogen atoms and adjacent slabs are separated by 
9.4 A of vacuum. Convergence required a kinetic-energy 
cutoff of 218 eV and a single (nonzero) k-point in the Bril- 
louin zone, yielding an accuracy of the calculated forces 
of 0.077 eV/A. In order to accelerate the convergence of 
the self-consistent cycle we introduced a metallic broad- 
ening of 0.054 eV. 

Figure [TJd) presents a normal incidence LEED pattern 
of the Si(331)-(12xl) reconstruction. The correspond- 
ing sketch containing reciprocal lattice vectors and spot 
labels {h^k) is shown in Fig. [It). The bulk-terminated 
surface is chosen as the reference for indexing. The re- 
ciprocal unit cell vectors Bi and B2 of the reconstructed 
surface can be expressed using the reciprocal unit cell 
vectors bi and b2 of the bulk-terminated (331) surface 
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However, since in general we find 11 satellite diffraction 
spots in between the integer spots along the [1 10] di- 
rection, the reconstruction is conventionally called the 
(12x1) reconstruction in the literature. 
The spot intensities in the LEED pattern in Fig. [TJd) ex- 
hibit a mirror symmetry along the [116] direction. Fur- 
thermore, a careful analysis of the LEED spot intensities 
as a function of energy (not shown) reveals that system- 
atically all half-order spots (±n/i/2,0), n being an odd 
integer, are missing for all beam energies (see empty cir- 
cles in Fig. [TJd) and c). Although spot intensities vary 
as a function of energy and even vanish at some energies 
due to diffraction effects, the (±n/i/2, 0) spots exhibit no 
intensity at all beam energies. Such missing spots in- 
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FIG. 2: (Color online) a) Large scale STM topography of the 
annealed Si(331) surface, b) High-resolution image with unit 
cell (full line) and the glide plane (dashed line) indicated. Bias 
voltage 2.0 V, set-point current 0.06 nA. 



dicate the existence of a glide plane along [116] in real 
space [13, [l8|. This glide plane also implies the mirror 
symmetry observed for the LEED intensities in reciprocal 
space [19]. Thus LEED firmly establishes the presence of 
a glide plane in the structure. 

We now turn to our STM data. The large scale to- 
pography image presented in Fig. [2^) shows the high 
quality of the surface. Besides the perfect long range 
order only a few local defects are present on the sur- 
face. As alread y o bserved by several other groups 
[ii,[20,[2l|,[23,[2i,[24], the Si(331)-(12xl) reconstruction 
consists of similar mounds arranged into zigzag chains 



running along the [116] direction separated by trenches. 
Focusing now on the high-resolution image in Fig. [2)3) we 
see that each of the mounds consists of five protrusions 
forming a pentagon. A further protrusion may be iden- 
tified linking two successive pentagons within the same 
chain. Pentagons with the same dimensions were al- 
ready observed onGe(110)-c(8xlO) [2o, z6] as well as on 
Ge(110)-(16x2) [25] and on Si(110)-(16x2) [9]. To our 
knowledge this is the first observation of such pentagons 
on a surface away from the (1 1 0) orientation^ indicating 
their high stability and confirming their fundamental role 
as an elementary building block in semiconductor surface 
reconstructions. 

Inspired by structural elements encountered on recon- 
structed Si(113) and Ge(113) surfaces [o, 7], An et al. 
[2.] have proposed an adatom-tetramer-interstitial (ATI) 
model for the pentagons observed on the (110) sur- 
faces. Its stability has subsequently been tested theo- 
retically by means of first-principles total-energy calcu- 
lations (Mini. 

In the following, we develop a coherent structural model 
for the Si(331)-(12x 1) reconstruction inspired by the ATI 
model and discuss similarities and differences with the 
(110) case. In a first step we need to determine the reg- 
istry of the surface reconstruction with respect to the 
bulk. This is not a trivial task, so we proceed in two 
steps. We focus first on the position of the surface recon- 
struction with respect to the underlying bulk along the 
[110] direction. Here the occurrence of the glide plane 
symmetry gives us the clue. Inspection of Fig. [2)3) shows 
that a glide plane is found at the center of the zigzag 
chain (dashed line) consistent with the observation of 
missing spots in the LEED pattern. The glide plane 
found on the surface must also be a glide plane of the 
bulk, since the space group of the bulk contains all sym- 
metry elements of the surface. A side view of the bulk 
terminated Si(331) surface is shown in Fig. [3^) with a 
top view of our model in Fig. 3b). The dashed line rep- 
resents the glide plane. Fig. [3t) offers a graphical proof 
for the existence of the glide plane in the bulk. After 
mirror refiection along the glide plane line, a translation 
by half the unit vector — Ai/2 is necessary to obtain the 
original registry. 

After determining the registry of the surface reconstruc- 
tion with respect to the [110] direction we now need to 
study the registry with respect to the [116] direction. 
Here we benefit from a comparison with the Si(llO) sur- 
face shown schematically in Fig. [3^) and f). For a sketch 
of the complete model for the Si(110)-(16x2) reconstruc- 
tion including the steps along [112] see Ref. [11]. Ac- 
cording to the ATI model for the Si(llO) surface, the 
pentagon seen in STM images consists of four adatoms 
(a,6,c,(i, empty circles) forming the tetramer and one sur- 
face atom (e, black dot) belonging to the first atomic 
layer (see Fig. [3]i)). The six- fold coordinated interstitial 
atom (/, empty circle) is located at the center of the pen- 
tagon slightly below the tetramer plane and consequently 
not directly visible in STM images (see simulated STM 
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FIG. 3: (Color online) a) Side and b) top view of the Si(331) 
surface. Black dots indicate the position of dangling bonds, 
some of which are saturated by pentamers and adatoms 
(empty circles) . The (1x1) (red) and (12x1) (blue) unit cells 
are shown. For comparison an experimental and a simulated 
STM image is underlaid. The dashed line represents the glide 
plane, c) Graphical proof for the existence of the glide plane. 

d) Sketch of the pentamer and the adatom with atom labels. 

e) Side and f) top view of the Si(llO) surface, (see text) 



images in Ref. [ul). The resulting structural element 
formed by the tetramer, atom e and / combined is called 
a pentamer [IT]. 
In order to integrate the pentamer building block into our 



model we note that the arrangement of danghng bonds on 
the bulk-truncated surface represented by black dots and 
marked by the double headed arrows in Fig. [3] differs be- 
tween the Si(llO) and Si(331) surface. Whereas dangling 
bonds on the Si(llO) surface occur in double rows run- 
ning along [110], double rows alternate with single rows 
of dangling bonds on the bulk-terminated Si(331) sur- 
face. For the Si(llO) structure, atom e actually belongs 
to one of these dangling bond double rows. Consequently 
we also anchor the two pentamers per unit cell required 
by STM on the double dangling bond row of the bulk- 
terminated Si(331) surface. It is important to note that 
anchoring the pentamers in this way provides exactly the 
same local binding configuration as on the Si(llO) sur- 
face, since careful comparison of the sideviews in Fig. 
[3ti) and e) reveals that the bulk-truncated Si(331) sur- 
face can be viewed as a highly stepped Si(llO) surface. 
The selected position of the pentagons agrees with the 
position observed in the STM images (see STM image 
behind the model in Fig. [3)3)). Note also that the sin- 
gle dangling bond row of the Si(331) surface is slightly 
lower-lying than the double dangling bond row, causing 
an inclination of the pentamer, in agreement with the 
lower intensity of the lobes associated with adatom b and 
c in the experimental STM image. 

Each pentamer saturates five of the surrounding dangling 
bonds [10]. By introducing two pentamers per Si(331)- 
(12x1) unit cell, the number of dangling bonds has been 
reduced from 36 to 26. Some of the remaining dangling 
bonds are saturated by simple adatoms as in the case of 
the Si(llO) surface. In the STM image in Fig. [2)d) the 
additional protrusion linking two successive pentamers 
indicates the location of a first adatom, labeled T4 in 
Fig. [3)3) in analogy with the convention used to label 
this adatom position on the (111) surface. The local 
binding configuration for this T4 adatom on the Si (331) 
surface is exactly the same as on the (111) surface, where 
adatoms are common structural building blocks. A fur- 
ther adatom labeled Hs saturates three more dangling 
bonds. 

Introducing two T4 and two Hs adatoms per unit cell 
into our structural model further reduces the number of 
dangling bonds from 26 to 14. For Si(llO) Stekolnikov 
et al. [ll| have noted that it is energetically more favor- 
able to leave some surface atoms, so called rest atoms, 
unsaturated than to introduce the maximum number of 
adatoms into the model, since this allows a further re- 
duction of the surface energy by electron transfer from 
the adatom to the rest atom in analogy with the Si(lll)- 
(7x7) case [13 • 

We have explicitly verified the stability of the reconstruc- 
tion containing two T4 and two H^ adatoms by means 
of first-principles calculations testing various alternative 
adatom configurations, among which the present one 
yields the highest stability (details on the theoretical in- 



vestigations will be published elsewhere [28]). In Fig. [3)3) 
we also show a simulated STM image based on relaxed 
coordinates of our structural model, obtained by inte- 
grating the local state density over a 2.0 eV energy win- 
dow corresponding to the experimental STM bias volt- 
age. The DFT image is in excellent agreement with the 
experimental STM image. The theoretical image repro- 
duces well the strong intensity of the lobes associated 
with atoms a, e and d and the weaker intensity associ- 
ated with atoms h and c of the pentamer. In general 
the topographical features in the simulated image show 
better resolution and less diffuse character than the cor- 
responding features observed in STM, since the simula- 
tion neglects any role of the tip structure in degrading 
image resolution [29]. In the experimental as well as in 
the simulated image, individual pentagons do not exhibit 
a mirror symmetry along the [116] direction indicating a 
distortion of the pentagon which allows to reduce internal 
strain. Furthermore the T4 and H^ adatoms in the simu- 
lated image are seen as a marked protrusion in agreement 
with experiment. 

It is interesting to note that on the Si(llO) surfaces, pen- 
tamers always occur in twin pairs rotated by 180^ with 
respect to each other, whereas on Si(331) all pentamers 
point into the same direction and are further apart from 
each other. Comparing the transition temperature of 
810^C [12] for Si(331)-(12xl) with 730^C [30] for Si(llO)- 
(16x2), it appears that Si (331) is the more favorable sur- 
face for pentamer formation (for comparison the tran- 
sition temperature of 870^C [3l| for the Si(lll)-(7x7) 
reconstruction is slightly higher). Further detailed total- 
energy calculations are required to quantify the trade- 
off between surface dangling bond reduction and induced 
surface stress and to compare the differences in the bond- 
ing configuration between the Si(331) and the Si(llO) 
surface. 

In summary, by combining STM, LEED, first-principles 
calculations and by comparing similarities and differ- 
ences between the Si(331)-(12xl) and Si(110)-(16x2) 
reconstructions, we have derived a complete structural 
model for the Si(331) surface containing the pentamer as 
an essential ingredient. Thus besides adatoms, dimers 
and tetramers, pentamers emerge as a universal building 
blocks for silicon surface reconstructions. 
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